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ABSTRACT 

Optimization  plays  a  key  role  in  environmental 
engineering  parameters  since  the  best  system 
performance  mainly  is  on  optimum  point  or  optimum 
range.  The  majority  of  wastewater  treatment  processes 
are  multi-variable  and  optimization  through  the 
classical  method  is  inflexible,  unreliable  and  time- 
consuming.  Thus,  response  surface  methodology 
(RSM),  as  a  very  efficient  design  and  widely  used 
technique,  can  be  adapted  for  optimization  of  various 
parameters  in  wastewater  treatment  processes.  RSM  is 
a  practical  mathematical  and  statistical  tool  that  can  be 
employed  for  analyzing  the  effects  of  several 
independent  factors  on  the  treatment  process  in  order  to 
obtain  the  maximum  benefit  from  the  process. 
Recently,  several  water  and  wastewater  treatment 
processes  have  been  optimized  for  treatment  different 
type  of  wastewaters  via  RSM  including;  textile  dye 
wastewater,  tannery  wastewater,  industrial  paint 
wastewater,  landfill  leachate,  olive  oil  wastewater,  and 
palm  oil  mill  effluent.  The  present  study  focuses  on  the 
usability  and  effectiveness  of  RSM  for  process 
parameters  modeling  and  optimization  in  wastewater 
treatment  studies.  In  this  paper,  some  of  the  RSM 
studies  were  reviewed  in  order  to  verify  the  usability  of 
RSM  and  its  limitations. 

Keywords:  Design  of  experiments,  Optimization,  Pulp 
and  Paper  Industry,  Response  Surface  Methodology, 
Wastewater. 


1.  Introduction 

The  paper  industry  is  one  of  the  largest  industries  in 
India,  consuming  large  amount  of  water  (Rajakumar  et 
al.  2011;  Bashir  et  al.  2012).  Indian  paper  mills  take  a 
wide  variety  of  cellulosic  and  non-cellulosic  raw 
materials.  The  pulp  and  paper  mill  in  India  utilizes 
different  cellulosic  materials  about  43%  wood  forest 
based,  28%  agro  based  and  29%  waste  paper  based  on 
the  total  installed  capacity  (Balakrishnan  1999).  These 
effluents  are  strongly  coloured  owing  to  the  presence  of 
lignin,  resin,  tannin  and  chlorophenolic  compounds  that 
are  resistant  to  biodegradation.  However,  very  few 
information  are  available  on  the  applicability  of  treating 
real  pulp  and  paper  mill  waste  water  by  anaerobic 
digestion  (especially  model/pilot/full-scale  HUASBR). 
Anaerobic  digestion  is  used  in  treatment  of  wastewater 
obtained  from  olive  mill  (Mohamed  1989),  dairy 
(Nemerow  et  al.  1991),  slaughter  house  (Mandal  and 
Bandana  1996)  and  dye  bath  effluents  (Carliell  et  al. 
1996).  But  in  case  of  Pulp  and  paper  mill  waste  water, 
it  is  not  used  as  widely  as  the  activated  sludge  process 
(Pokhrel  and  Viraraghavan  2004).  Hence,  pulp  and 
paper  mill  waste  water  needs  detailed  study.  The  Up 
flow  anaerobic  sludge  blanket  (UASB)  reactor 
introduced  by  (Lettinga  et  al.  1980)  has  become  a 
popular  high-rate  anaerobic  treatment  system 
throughout  the  world  (Lettinga  et  al.  1993).  In  the 
journal  article  (Rindone  et  al.  1991),  the  Hybrid  up 
flow  anaerobic  sludge  blanket  (HUASB)  reactor  is  a 
new  concept  which  is  the  hybridized  version  of  an 
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UASB  reactor  with  a  random  packing  media  at  the  top 
of  the  reactor.  In  which  upper  portion  of  50-70%  is 
filled  with  either  floating  (or)  stationary  materials  to 
retain  some  of  the  escaping  biomass.  This  HUASB 
reactor  combines  the  merits  of  UASB  (Suspended 
growth  system)  and  fixed  film  reactors  (Attached 
growth  system).  This  type  of  reactor  is  of  particular 
value  in  a  situation  when  the  rate  of  sludge  granulation 
is  slow  and  there  is  a  need  to  accelerate  the  reactor 
start-up  (Badroldin  et  al.  2008).  In  recent  years, 
HUASB  reactor  has  proved  to  be  a  more  versatile 
anaerobic  treatment  in  treatment  of  industrial 
wastewater.  Therefore  an  attempt  will  be  taken  to 
measure  the  treatability  performance  of  HUASB  reactor 
using  pulp  and  paper  mill  wastewater  (Hemalatha, 
Sanchitha  and  Keerthinarayana,  2014). 

2.  Microbial  Activities  in  the  Anaerobic  Process 

The  number  and  types  of  microorganisms  present  in 
digesters  are  likely  to  depend  upon  the  type  of  digester, 
its  operating  conditions  and  the  waste  composition.  The 
metabolic  stages  involved  in  the  production  of  methane 
from  waste  in  Anaerobic  Digestion  occur  in  4  distinct 
processes  (Fig.  1); 

2.1.  Hydrolysis 

Complex  organic  matter  is  decomposed  into  simple 
soluble  organic  molecules  using  water  to  split  the 
chemical  bonds  between  the  substances.  This  is  where 
solid  complex  organics,  cellulose  proteins,  lignins,  and 
lipids  are  broken  down  into  soluble  (liquid)  organic 
fatty  acids.  The  results  are  soluble  monomers. 
Hydrolytic  bacteria  are  responsible  for  the  creation  of 
monomers.  Enzymes  excreted  from  the  bacteria,  such 
as  cellulase,  protease,  and  lipase,  catalyse  hydrolysis. 
Therefore,  the  more  complex  the  feedstock  then  the 
hydrolytic  phase  is  relatively  slow.  A  hydrolysis 
reaction  where  organic  waste  is  broken  down  into  a 
simple  sugar,  in  this  case  glucose  can  be  seen  in  the 
following  Eq.  1 . 

C6H10O4+  2H20  C6Hi206  +  2H2 

2.2.  Fermentation/Acidogenesis 

The  chemical  decomposition  of  Carbohydrates,  proteins 
and  fats  by  enzymes,  bacteria,  yeasts  in  the  absence  of 
oxygen.  Hydrolysis  is  immediately  followed  by  the 
acid- forming  phase  of  Acidogenesis.  Here  acidogenic 
bacteria  turn  products  of  hydrolysis  into  mostly  short 
chain  (volatile)  acids  (e.g.  formic  or  lactic),  ketones 
(e.g.  ethanol  or  acetone)  and  alcohols.  The  specific 
concentrations  of  products  formed  here  vary  with  the 


type  of  bacteria,  culture  conditions,  such  as  temperature 
and  pH. 

2.3.  Acetogenesis 

The  fermentation  products  are  converted  into  acetate 
hydrogen  and  carbon  dioxide  by  so  called  acetogenic 
bacteria.  Here  the  Biological  Oxygen  Demand  (BOD) 
and  Chemical  Oxygen  Demand  (COD)  are  reduced. 
Acetogenesis  occurs  through  carbohydrate  fermentation 
which  acetate  is  the  main  product,  and  other  metabolic 
processes.  In  Eq.2  reaction  in  the  acid- forming  stages  is 
shown  below,  glucose  is  converted  to  ethanol. 

C6Hi206  ~  2CH3CH2OH  +2C02 

2.4.  Methanogenesis 

Methane  (CH4)  is  formed  from  acetate  and 
hydrogen/carbon  dioxide  by  methanogenic  bacteria. 
Most  methanogenic  bacteria  utilize  H2  &  C02,  but 
species  of  only  two  genera  Methanosarcina  and 
Methanorthrix,  can  produce  methane  from  acetic  acid. 
The  acetogenic  bacteria  grow  in  close  association  with 
the  methanogenic  bacteria  during  the  4th  stage  process 
(Osman,  2014). 

2.5  Advantages  of  Anaerobic  fermentation 

1.  It  can  treat  a  wide  range  of  organic  wastes 
including  industrial  wastewater 

2.  No  odour  nuisance  during  the  process  and  the 
reduction  approximately  80%  of  the  odour 
potential. 

3.  It  is  relatively  small  in  size  to  the  amount  of 
waste  treated. 

4.  Anaerobic  fermentation  projects  can  directly 
boost  the  local  rural  economy  through  creating 
jobs  in  the  anaerobic  fermentation  development 
and  indirectly  through  increasing  disposable 
income  in  rural  areas. 

5.  Reducing  land  and  waste  pollution:  Poor 
disposal  of  wastewater  can  cause  land  and 
ground  water  pollution.  Anaerobic  fermentation 
creates  integrated  management  system  which 
reduces  the  likelihood  of  this  happening  and 
reduces  the  likelihood  of  fines  been  imposed  of 
such  pollution. 

6.  Anaerobic  fermentation  provides  onsite  energy 
for  the  process 

The  biogas  produced  is  a  renewable  energy  source  and 
used  as  a  transport  fuel  or  to  produce  electricity 
displaces  fossil  fuel  energy  and  they’re  by  reducing  the 
emissions  of  green  house  and  acidifying  gasses 
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1.  Anaerobic  digestion  is  also  a  technology  that  can  make  a  significant  contribution  to  the  management  of 
organic  waste. 

2.  The  fiber  form  the  anaerobic  digestion  process  can  be  used  as  a  good  soil  conditioner  and  the  liquor  can  be 
used  as  a  fertilizer 


Figure  1 :  Process  flow  of  anaerobic  digestion  mechanism 


3.  Upflow  anaerobic  sludge  blanket  reactor 
(UASBR) 

One  of  the  most  notable  developments  in  anaerobic 
treatment  process  technology  is  the  upflow  anaerobic 
sludge  blanket  (UASB)  reactor  (Lettinga  et  al.,  1980) 
with  its  wide  applications  in  relatively  dilute  municipal 
wastewater  treatment  (Tchobanoglous  et  al.,  2003). 
Influent  flow  distributed  at  the  bottom  of  the  UASB 
reactor  travels  in  an  upflow  mode  through  the  sludge 
blanket  and  passes  out  around  the  edges  of  a  funnel 
which  provides  a  greater  area  for  the  effluent  with  the 
reduction  in  the  upflow  velocity,  enhancement  in  the 
solids  retention  in  the  reactor  and  efficiency  in  the 
solids  separation  from  the  outward  flowing  wastewater 
(Ersahin  and  Ozgun,  no  date).  Granules  which  naturally 
form  after  several  weeks  of  the  reactor  operation 
consist  primarily  of  a  dense  mixed  population  of 
bacteria  that  is  responsible  for  the  overall  methane 
fermentation  of  substrates  (Rittmann  &  McCarty, 
2001).  Good  settleability,  low  retention  times, 
elimination  of  the  packing  material  cost,  high  biomass 
concentrations  (30000-80000  mg/L),  excellent 
solids/liquid  separation  and  operation  at  very  high 
loading  rates  can  be  achieved  by  UASB  systems 
(Speece,  1996).  The  only  limitation  of  this  process  is 
related  to  the  wastewaters  having  high  solid  content 


which  prevents  the  dense  granular  sludge  development 
(Tchobanoglous  et  al.,  2003).  Design  OLR  is  typically 
in  the  range  of  4  to  15  kg  COD/m3.day  (Rittmann  & 
McCarty,  2001). 

4.  Major  pulp-and-paper  processes 

The  main  processes  in  pulp  manufacturing  and  paper 
making  are  divided  into  five  major  groups  including 
mechanical,  chemical,  chemo-mechanical,  and  thermo¬ 
mechanical  pulping  as  well  as  papermaking,  each 
producing  a  high  volume  of  wastewater  with  special 
characteristics. 

4.1.  Wastewater  production  in  the  pulp-and-paper 
processes 

Wood  preparation,  pulp  washing,  pulp  bleaching  and 
paper  making  processes  as  well  as  the  digester  house 
are  the  major  wastewater  producers  in  the  pulp-and- 
paper  industry  (The  World  Bank  Group,  1999).  The 
generated  wastewaters  have  a  high  content  of  BOD  and 
various  concentrations  of  other  contaminants  that 
depend  on  the  types  of  applied  processes.  For  instance, 
wood  preparation  wastewater  has  suspended  solids, 
BOD,  dirt,  and  fibers  while  the  produced  wastewater  in 
the  digesters  house  contains  resins,  fatty  acids,  color, 
BOD,  COD,  AOX,  and  VOCs  (Pokhrel  and 
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Viraraghavan,  2004).  The  magnitude  of  GHG  emission 
in  the  WWTPs  depend  on  the  concentrations  of 
nutrients  in  wastewater,  reactor  operating  temperature, 
type  of  treatment  processes  employed  and  the 
contaminant  removal  efficiency(Ashrafi  et  al.,  2013b). 

4.2.  Wastewater  treatment 

The  growing  concerns  about  the  use  of  fresh  water, 
increased  economic  considerations  and  stringent 
environmental  regulations  have  highlighted  the 
importance  of  water  use  and  efficient  wastewater 
treatment  in  the  pulp-and-paper  industry.  The  recycling 
and  reuse  of  the  generated  wastewater  after  its  proper 
treatment  is  the  key  to  reduce  fresh  water  use.  This 
procedure  will  minimize  external  discharges  to  the 
environment  while  advancing  environmental 
conservation  by  reducing  fresh  water  consumption.  In 
view  of  this,  the  main  objective  of  WWTPs  is  to 
remove  contaminants  from  wastewater  using  a  series  of 
physicochemical,  biological,  and  integrated  treatment 
processes.  Ashrafi  et  al.  (2013a;  2013b)  developed 
comprehensive  mathematical  models  (steady-state  and 
dynamic  models)  to  estimate  the  generated  GHG 
emission  by  WWTPs.  In  the  development  of  these 
models,  it  was  assumed  that  the  bioreactors  are 
completely  mixed  (Ashra,  Yerushalmi  and  Haghighat, 
2015): 

Accumulation  =  influent  -  effluent  +  net  growth 
(production  by  reaction) 

5.  Response  Surface  Methodology 

The  Response  Surface  Methodology  (RSM)  emerged  in 
the  1950s  (BOX  et  al.,  1951;  BOX 
et  al.,  1951)  within  the  context  of  Environmental 
Engineering  in  an  attempt  to  construct  empirical 
models  able  to  find  useful  statistical  relationships 
between  all  the  variables  making  up  an  industrial 
system.  Optimization  plays  a  key  role  in  this  domain 
and  parameters  to  optimize  the  system’s  performance. 
Majority  of  the  wastewater  treatment  processes  are 
multi-variable  and  optimization  through  the  classical 
method  is  inflexible,  unreliable  and  time-consuming. 
Thus,  an  alternative  method  which  will  be  more 
effective  and  can  be  adapted  for  parameter  optimization 
of  various  wastewater  treatment  processes  is  favored. 
Response  surface  methodology  (RSM)  is  one  of  the 
most  efficient  and  widely  used  mathematical  and 
statistical  tools  for  system  performance  optimization. 
RSM  can  be  employed  to  optimize  and  analyze  the 
effects  of  several  independent  factors  on  a  treatment 
process  to  obtain  the  maximum  output.  This 


methodology  is  based  on  experimental  design  with  the 
final  goal  of  evaluating  optimal  functioning  of 
industrial  facilities,  using  minimum  experimental 
effort.  Here,  the  inputs  are  called  factors  or  variables 
and  the  outputs  represent  the  response  that  generates 
the  system  under  the  causal  action  of  the  factors. 
Afterwards,  the  use  of  the  RSM  was  shown  in  the 
design  of  new  processes  and  products.  In  recent  years  it 
is  being  applied  successfully  in  other  scientific  fields 
such  as  biology,  medicine,  and  economy  (Jimenez- 
contreras  and  Torres-salinas,  2009).  MYERS  et  al. 
(2004)  has  exhaustively  reviewed  the  literature  in  the 
sense,  describing  the  developments  and  applications  of 
this  methodology.  Very  recently,  RSM  has  been  used 
even  to  validate  new  experimental  methods  (JURADO 
et  al.,  2003). 

Besides  analyzing  the  independent  variables  effects, 
this  experimental  methodology  also  generates  a 
mathematical  model.  The  graphical  viewpoint  of  the 
mathematical  model  has  led  to  the  term  RSM.  The 
relationship  between  the  responses  and  the  inputs  is 
given  in  Eq.  (3): 

Y=f(XI,X2,X3,....,Xn)±s 

Where  Y  is  the  response,  f  is  the  unknown  function  of 

response,  xi,  x2,  . xn  are  the  input  variables  which 

can  affect  the  response,  n  is  the  number  of  the 
independent  variables  and  e  is  the  statistical  error  that 
represents  other  sources  of  variability  not  accounted  for 
by  f.  After  selection  of  the  design,  the  model  equation 
is  defined  and  coeffecnts  of  the  model  equation  are 
predicted. 

As  the  result  of  sequential  model  sum  of  squares 
suggestion  of  software  and  software  suggestion,  the 
quadratic  model  was  selected.  In  the  case  that  total 
number  of  experiments  is  n;  the  response  surface  can  be 
expressed  as  follows  using  matrix  notation  of  the 
model. 

Y  =  Xp  ±  e 

Where 
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Where  s  is  random  error. 

The  following  equation  is  for  calculating  the  total 
number  of  experiments 

N  =  n3  +  2n  +  nc 

Where  N  is  the  total  number  of  experiments  and  n  is  the 
number  of  factors. 

The  selected  independent  variables  were  coded 
according  to  equation  (4) 

X  —  X 

X.  =  *  0  i  - 1, 2, . k 

AX 

Where  X;  refers  to  coded  value  of  the  ith  independent 
variable,  X0  is  the  value  of  X,  at  the  center  point  and 
AX  is  the  step  change  value  [12]  the  response  function 
(Y)  is  measured  as  the  yield  in  processes.  The  model 
used  in  RSM  is  commonly  a  full  quadratic  equation  in 
the  diminished  form  of  this  equation.  The  second  order 
model  can  be  written  as  follows. 


Y  =  fio  +  Z  fix  +  Z  fix  +  Z  fiiXXij  +  e 

i=l  i=l  i*  y=l 

Where  (30  is  the  value  of  the  fixed  response  at  the 
center  point  of  the  design;  (3i,  (3ii,  and  (3ij  are  the  linear, 
quadratic  and  interaction  effect  regression  terms, 
respectively;  xi  denotes  the  level  of  the  independent 
variable;  n  is  the  number  of  independent  variables;  and 
s  is  random  error.  Normally,  a  second  order  polynomial 
equation  was  derived  as: 

Y  —  /?(l  +  +  fiXi  fix 3  fi\XXi  fiiXXi  fix  i 

5.1.  RSM  application  in  optimization  of  wastewater 
treatment  processes 

In  this  section,  several  published  RSM  studies  in  the 
last  few  years  were  reviewed  focusing  on  the  usability 
of  RSM  for  optimization  of  various  types  of  wastewater 
treatment  processes.  Table  1  summarizes  different 
applications  of  RSM  in  optimizing  various  types  of 
wastewater  treatment  techniques. 


Table  1:  Various  application  of  RSM  in  wastewater  treatment  process  optimization 


Wastewater 

type 

Treatment 

technique 

Independent  variables 

Responses 

Ref. 

Textile 

wastewater 

Electrochemical 

Oxidation 

pollution  load, 
applied  potential,  electrolyte 
concentration, 

temperature,  and  reaction  time 

COD,  color  turbidity 
removals 

[19] 

Textile  dyes 

Adsorption 

Initial  dye  concentration,  initial 
solution  pH  and 
temperature 

the  amount  of  dye 
adsorbed  at 
equilibrium 

[39] 

Landfill 

leachate 

Coagulation- 

Flocculation 

coagulant  dosage  and  pH 

COD,  turbidity,  color, 
and  TSS  removals 

[15] 

Landfill 

leachate 

Fenton  Oxidation 

pH,  Reaction  time,  Initial 
concentrations  of  H202,  and  Ferrous 
ion  concentration 

COD,  color  and 
iron  removals 

[25] 

Landfill 

Leachate 

Electrochemical 

oxidation 

electrolyte 

concentrations,  current  density  and 
reaction  time 

COD  and  color 
removals 

[26] 

Landfill 

leachate 

Ion  Exchange 
(cation  resin) 

Cation  resin  dosage,  contact  time, 
shaking  speed 

Ammonia  removal 

[9] 

Landfill 

leachate 

Ion  Exchange 
(anion  resin) 

anion  resin  dosage,  contact  time, 
shaking  speed,  pH 

COD,  color,  turbidity, 
and  SS  removals 

[10] 

Pulp  mill 

wastewater 

Coagulation- 

Flocculation 

coagulant  dosage,  flocculant  dosage 
and  pH 

turbidity  removal, 
lignin  removal  and 
clean  water  recovery 

[47] 

Slaughterhouse 

Wastewater 

Electrochemical 

oxidation 

current  density,  reaction  time  and 
influent  COD 

COD,  BOD,  and  color, 
removal  efficiencies 
and  effluent  pH 

[5] 
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Simulated 

industrial 

wastewater 

Photo-Fenton 

process 

initial  pH  values,  concentration  of 
iron  catalyst 

and  the  concentration  of  H202  type  of 
UV  irradiation 

Mineralization  rate 

[13] 

Dairy 

wastewater 

Electrochemical 

oxidation 

current  density,  dosage  of  sodium 
chloride  (NaCl),  electrolysis  time  and 
pH, 

COD  removal 

[20] 

Petroleum 

refinery 

effluent 

Upflow  anaerobic 
sludge  blanket 
(UASB)  bioreactor 

hydraulic  retention  time  (HRT) 
influent  COD,  up  flow  velocity  (Vup) 

COD  removal,  and 
rate  of  biogas 
production 

[35] 

Oily 

wastewater 

emulsion 

Electro¬ 

coagulation 

Current  density,  pH,  and 
Electrocoagulation  time 

turbidity  and  COD 
removals 

[45] 

Drinking  water 

Coagulation- 

flocculation 

coagulant  dose  and  pH 

removal  efficiency  of 
turbidity 

and  dissolved  organic 
carbon  (DOC) 

[46] 

Olive  oil  mill 
wastewater 

Fenton’s 

peroxidation 

he  ratio  of  hydrogen  peroxide-to- 
Fe(II) .  Fe(II)  concentration  ,  and 
H202  concentration 

COD,  total  phenolics 
(TP),  color  and 
aromatocity  removal 

[1] 

Table  2:  Feasibility  of  different  treatment  method  for  different  wastewater  sources 


Wastewater  type 

Independent  variables 

Responses 

Reactor 

Start  up 
period 

Typical 
loading 
rates  (kg 
COD/m3 
day) 

HRT 

(d) 

Textile  wastewater 

Pollution  load, 
temperature,  and  reaction 
time 

COD,  color 
turbidity 
removals 

CSTR 

0.25-3 

10-60 

Pulp  mill  wastewater 

Pollution  load, 
temperature,  and  reaction 
time 

turbidity 

removal, 

lignin 

removal  and 
clean  water 
recovery 

UASBR 

4-16 

10-30 

0.5-7 

Dairy  wastewater 

current  density,  dosage  of 
sodium  chloride  (NaCl), 
electrolysis  time  and  pH 

COD 

removal 

AFB 

3-4 

1-100 

0.2-5 

Subjected  to  field 
conditions 

6.  CONCLUSION 

Response  surface  methodology  (RSM)  is  a  widely  used 
technique  which  can  be  adapted  for  parameters 
optimization.  Statistical  and  diagnostics  analysis 
indicated  that  RSM  is  a  reliable  tool  to  optimize  of 


process  parameters.  RSM  is  successfully  applied  for 
optimization  of  several  wastewater  treatment  processes. 
In  the  current  study,  s  RSM  method  and  some  of  its 
applications  that  published  recently  were  reviewed  in 
order  to  verify  the  usability  of  RSM  and  its  limitations. 
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It  was  observed  that  RSM  has  several  advantages  in 
wastewater  treatment  process  optimization.  However, 
in  most  of  the  reviewed  articles  that  implement  RSM 
for  optimization  of  wastewater  treatment  processes, 
there  is  not  sufficient  or  adequate  preliminary  work  in 
regards  to  the  range  of  process  variables.  This  could 
lead  to  inaccurate  outcome. 
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